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Путем анализа экспериментальных данных исследуется коллективная динамика ан- 
самблей из нескольких десятков нефронов, расположенных на поверхности почки. На 
основе вейвлет-анализа изучается эффект захвата мгновенных частот и фаз колебаний, 
обусловленных механизмом канальцево-гломерулярной обратной связи. Показано, что в 
формировании кластеров синхронной динамики участвуют структурные элементы поч- 
ки, относящиеся к разным нефронным «деревьям», а подстройка частот и фаз колебаний 
больших групп нефронов наблюдается только на определенных фрагментах эксперимен- 
тальных записей. Установлено, что значительные группы нефронов, расположенных на 
разных участках поверхности почки, демонстрируют эффект синфазной синхронизации 
колебаний. 
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Введение 

Синхронизация автоколебаний относится к числу фундаментальных явлений 
природы [1-6], многочисленные примеры которых хорошо известны в динамике 
биологических систем. В их числе можно упомянуть синхронизацию электрической 
активности нейронных ансамблей [11], синхронизацию сокращений сердца и рит- 
ма дыхания [21], вынужденную синхронизацию сердечного ритма периодическим 
звуковым или световым воздействием [22], кооперативную динамику популяций на- 
секомых [19] и т.д. 

Одним из примеров проявления эффекта синхронизации в функционировании 
объектов живой природы служит динамика нефронов - структурных элементов поч- 
ки. Нефрон представляет собой объект размером порядка 100 мкм, принимающий 
участие в фильтрации крови и регуляции кровяного давления [10]. Для сравнения 
можно отметить, что почка человека содержит примерно миллион, а почка крысы 
- около тридцати тысяч нефронов. Расположение нефронов визуально напоминает 
дерево, на «ветках» которого обычно находятся от одного до трех таких элементов. 
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Организация нефронных «деревьев» в виде параллельных структур приводит к то- 
му, что общий почечный кровоток делится между большим числом нефронов, и на 
каждый из них приходится только очень малая часть потока крови, поступающего 
через почечную артерию [11, 12]. 

Известно несколько механизмов авторегуляции почечного кровотока на уровне 
отдельных нефронов, среди которых выделяется канальцево-гломерулярная обрат- 
ная связь (КГОС) [13, 14]. Данный механизм приводит к генерации автоколебаний 
давления жидкости с периодом примерно 30-40 секунд, которые регистрируются 
как в малых кровеносных сосудах (артериолах), так и в фильтрате, протекающем 
по канальцам нефрона. Согласно результатам экспериментальных исследований на 
крысах, соответствующие колебания различаются в норме и при гипертонии - при 
нормальном артериальном давлении они являются почти периодическими, а при по- 
вышенном давлении (генетическая форма гипертонии или искусственно вызванная 
гипертония Голдблетта) эти колебания становятся сильно нерегулярными (хаотиче- 
скими) [15-18]. В связи с этим изучение динамики нефронов представляет собой 
одну из прикладных задач нелинейной динамики, а выявление причин хаотизации 
колебаний вызывает несомненный интерес для более глубокого понимания механиз- 
мов, приводящих к развитию почечной гипертонии. 

Недавно проводившиеся исследования [19-21] продемонстрировали, что эф- 
фекты синхронизации колебаний давления жидкости в канальцах соседних структур- 
ных элементов почки в норме и при гипертонии различаются. Так, в случае патоло- 
гии длительность участков захвата частот или фаз колебаний уменьшается примерно 
в 3-4 раза [21], а полная синхронизация колебаний парных нефронов сменяется ча- 
стичной синхронизацией [19,20]. Однако до настоящего времени взаимная динамика 
изучалась лишь для случая малых ансамблей - парных нефронов или триплетов. Это 
было связано с двумя причинами - сложностью проведения экспериментов по одно- 
временной регистрации сигналов большого числа структурных элементов почки и 
теоретическими представлениями, согласно которым синхронизация должна наблю- 
даться для нефронов, относящихся к одной междольковой артерии. Коллективная 
динамика больших групп нефронов до сих пор не рассматривалась, поэтому су- 
ществует целый ряд открытых вопросов о взаимодействии структурных элементов 
почки, относящихся к разным нефронным «деревьям». Данная работа представляет 
собой одну из первых попыток изучения синхронизации ансамблей из нескольких 
десятков корковых нефронов, регистрируемых на поверхности почки крысы. Заме- 
тим, что проведение экспериментов іп ѵіѵо возможно только на корковых нефронах, 
однако существующие сейчас экспериментальные возможности позволяют рассмат- 
ривать динамику структурных элементов, расположенных на значительной части по- 
верхности почки. 


1. Методика исследования 

1.1. Регистрация сигналов нефронных ансамблей. Для изучения коллек- 
тивного поведения больших ансамблей нефронов в Институте Панум (Университет 
Копенгагена) были проведены несколько серий экспериментов на крысах, в рамках 
которых осуществлялась регистрация динамики структурных элементов на поверх- 


15 



ности почки. Первоначально регистрация проводилась с помощью инфракрасной ка- 
меры, позволяющей по изменению температурного режима анализировать скорость 
потока крови в артериолах. Однако применение такого метода не обеспечивало вы- 
сокого качества регистрации динамики кровотока в микрососудах. Кроме того, по 
полученной видеозаписи температурного режима нельзя было установить, где на- 
ходятся нефроны. Чтобы идентифицировать расположение структурных элементов 
на поверхности почки, с использованием микроскопа делалась фотография поверх- 
ности, и путем наложения фотографии на изображение, полученное с видеокамеры, 
можно было установить координаты расположения нефронов (в пикселях). Тем са- 
мым получалась информация о том, какие координаты выбрать для извлечения из ви- 
деофайла временной динамики, соответствующей каждому обнаруженному нефрону. 
Тем не менее, в связи с невысокой точностью данного способа регистрации коллек- 
тивной динамики структурных элементов почки, в дальнейшем он не использовался. 

В рамках последующих экспериментов применялся другой метод - лазерной 
спекл-интерферометрии (ЛСИ, в зарубежной литературе Іаяег .чрескіе іта§іп§), кото- 
рый обеспечивает высокое пространственно-временное разрешение при регистрации 
скорости микрокровотока в распределенных в некоторой области пространства кро- 
веносных сосудах [22,23]. Этот метод имеет преимущество по сравнению с часто 
используемым в физиологических исследованиях методом лазерной доплеровской 
флуометрии (ЛДФ) [24]. Последний подход требует проведения сканирования для 
осуществления пространственного разрешения, а механическая процедура сканиро- 
вания поверхности, представляющей интерес для исследования, ограничивает разре- 
шающие способности метода ЛДФ. Сходные проблемы отмечались при рассмотре- 
нии альтернативного способа - изучения спеклов, меняющихся во времени [25,26]. 
Метод лазерной спекл-интерферометрии к настоящему времени продемонстрировал 
свою эффективность при изучении кровотока в сетчатке глаза [27], мониторинге це- 
ребрального кровотока [28] и т.д. По аналогии с использованием инфракрасной ка- 
меры, данные лазерной спекл-интерферометрии дополнялись фотографией верхней 
поверхности почки, позволяющей установить расположение нефронов. 

Экспериментальная база данных, записанных в институте Панум с применени- 
ем ЛСИ-метода, включала сигналы ансамблей нефронов пятнадцати крыс. На верх- 
ней поверхности почки каждой крысы идентифицировалось от 40 до 100 нефронов в 
зависимости от эксперимента (для анализа выбирались только данные четко разли- 
чимых нефронов). Временное разрешение метода составляло 1 Гц, что обеспечивало 
получение порядка 30-40 отсчетов на один характерный период колебаний, обуслов- 
ленных механизмом КГОС. Длительность каждой записи составляла около 30 минут, 
в результате для каждого нефрона регистрировалось до 100 периодов ритма КГОС. 

1.2. Анализ экспериментальных данных. Записанные экспериментальные 
данные анализировались на основе непрерывного вейвлет-преобразования, которое 
обеспечивает возможность изучения локальных изменений структуры сигналов [29, 
30]. Если классический спектральный анализ дает информацию о частотном соста- 
ве рассматриваемого процесса и характерных ритмах колебаний, то вейвлет-анализ 
позволяет извлекать информацию о таких характеристиках, как мгновенные ампли- 
туды, мгновенные частоты и фазы ритмических процессов. Особенностями вейвлет- 
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анализа является использование базисных функций, локализованных во временной и 
в частотной областях, и наличие подвижного частотно-временного окна, меняюще- 
гося с масштабом наблюдения. Последнее отличает вейвлет-преобразование от окон- 
ного преобразования Фурье - вейвлетная функция имеет разную протяженность на 
низких и высоких частотах, и в этом случае быстрые изменения структуры сигнала 
анализируются по коротким фрагментам экспериментальной записи, а медленные - 
по более длительным. 

Сравнительному исследованию фурье- и вейвлет-преобразований, а также де- 
талям процедуры вейвлет-анализа сигналов посвящено большое количество моно- 
графий (например, [29-32]). Более того, вейвлеты стали восприниматься как стан- 
дартный метод цифровой обработки сигналов в самых разных областях науки и 
техники (свыше 3 миллионов ссылок в сети Интернет), поэтому в данной статье 
не будем останавливаться на деталях вейвлет-преобразования. Отметим лишь, что в 
качестве базисной функции в работе использовался вейвлет Морле 

т = ^ е~ і2 !\ 

центральная частота которого /о была выбрана равной единице. Процедура непре- 
рывного вейвлет-преобразования и идентификации мгновенных частот и амплитуд 
колебательных процессов детально рассмотрена в [29, 30], она неоднократно приме- 
нялась и обсуждалась в наших предыдущих работах (например, [33,34]). 

2. Результаты 

Прежде чем приступать к рассмотрению кооперативного поведения больших 
ансамблей нефронов, обратимся к случаю динамики нескольких соседних структур- 
ных элементов почки. Как показано в работах [20,21], для динамики взаимодей- 
ствующих нефронов нормотензивных крыс характерны участки синхронизации, ко- 
торые включают в среднем примерно 10-12 периодов ритма КГОС, причем фрагмен- 
ты синхронных колебаний сменяются 
фрагментами несинхронного режима. 

Результаты отмеченных работ [20, 21] 
были получены по экспериментальным 
данным зависимостей давления в прок- 
симальных канальцах нефронов (кате- 
терный метод регистрации, обеспечи- 
вающий сравнительно низкий уровень 
помех). Чтобы продемонстрировать, что 
диагностика эффекта синхронизации не 
зависит от способа регистрации экспе- 
риментальных данных, на рис. 1 пред- 
ставлены временные зависимости мгно- 
венных частот колебаний, вычисленные 



Рис. 1. Динамика мгновенных частот колебаний, 
обусловленных механизмом КГОС, для трех нефро- 
нов, расположенных на поверхности почки 
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по экспериментальным записям, полученным в рамках ЛСИ-метода. Для нагляд- 
ности на этом рисунке изображены временные зависимости мгновенных частот 
только трех нефронов из числа тех, для которых наблюдаются синхронные ко- 
лебания. Первоначально (в диапазоне примерно до 600 секунд) регистрируются 
случайные изменения мгновенных частот. Вторая часть представленных зависимо- 
стей (в диапазоне от 600 до 1200 секунд) соответствует случаю синхронной ди- 
намики, когда мгновенные частоты колебаний всех трех нефронов оказываются 
захваченными. 

Наличие синхронизации колебаний малых групп соседних нефронов является 
ожидаемым результатом, более интересен ответ на другой вопрос - сколько нефро- 
нов могут синхронизовать свои колебания. Проводя изучение динамики структур- 
ных элементов на поверхности почки, мы лишены возможности заглянуть внутрь 
и изучать поведение нефронов, удаленных от поверхности. Мы можем лишь пред- 
полагать, что эффект синхронизации должен наблюдаться для нефронов в преде- 
лах одного «дерева», и эти предположения базируются как на общих соображениях 
(элементы, относящиеся к одному кровеносному сосуду, должны взаимодействовать 
между собой), так и на данных компьютерного моделирования [12]. Однако наличие 
большого числа корковых нефронов обеспечивает возможность изучения количества 
кластеров, которые можно идентифицировать на поверхности почки, а также опре- 
деления числа структурных элементов, участвующих в формировании кластеров. 

Проводя изучение динамики нефронных ансамблей, мы получили довольно 
неожиданный результат, состоящий в том, что захват частот колебаний наблюдается 
не только для соседних нефронов, но и для структурных элементов почки, находя- 
щихся на значительном расстоянии друг от друга и относящихся к разным нефрон- 
ным деревьям! Один из примеров представлен на рис. 2, где разными символами 
отмечены элементы, принадлежащие двум разным кластерам синхронизации. В дан- 
ном случае выбранное разрешение (600x600 пикселей) соответствует примерно 80% 
верхней поверхности почки. 

Как видно из приведенного ри- 
сунка, кластеры синхронной динами- 
ки нефронов не являются локализован- 
ными на поверхности почки и могут 
включать структурные элементы, рас- 
положенные на значительном расстоя- 
нии друг от друга. Кроме того, дан- 
ные кластеры включают большое коли- 
чество нефронов, а не малые ансамбли 
из нескольких элементов. Пример, при- 
веденный на рис. 2, является нагляд- 
ным, но не очень типичным - лишь 
в трех экспериментах из пятнадцати 

Рис. 2. Пример идентификации двух кластеров син- 
хронной динамики, соответствующих колебаниям с было зафиксировано наличие двух СО- 
частотами примерно 0.017 и 0.019 Гц существующих синхронных кластеров. 
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Более типичной является ситуация, ко- 
гда идентифицируется только один кла- 
стер и часть нефронов, которые не де- 
монстрируют синхронной динамики. В 
ходе проведенного статистического ана- 
лиза нами было установлено, что от 10 
до 60% корковых нефронов принима- 
ет участие в формировании синхронно- 
го кластера (в среднем около 40%). При 
этом кластер может существовать срав- 
нительно небольшое время (около 8- 
10 периодов колебаний), после чего его 
размер уменьшается и часть нефронов 
выходит из синхронного режима. 



Рис. 3. Длительности участков синхронизации для 
каждого индивидуального нефрона (черные линии). 
По оси ординат отложен номер нефрона 


Так как коллективная динамика структурных элементов почки носит нестаци- 
онарный характер, наряду с усредненными статистическими данными необходимо 
изучать временную динамику ансамбля. Это можно сделать различными способа- 
ми, и один из возможных вариантов приведен на рис. 3. На рисунке изображена 
динамика 40 выделенных нефронов с поверхности почки. Черными линиями пока- 
зан случай синхронного режима, а именно ситуация, когда каждый индивидуальный 
нефрон синхронизирует колебания по крайней мере еще с одним нефроном в течение 
не менее 10 периодов колебаний (при этом возможно, что синхронизация наблюда- 
ется с несколькими нефронами или даже несколькими десятками нефронов). Как 
следует из рисунка, в разные моменты времени в формировании синхронного кла- 
стера могут участвовать разные нефроны. Это позволяет выдвинуть гипотезу о том, 
что синхронизация нефронных ансамблей обусловлена не только геометрией «де- 
ревьев», а скорее, представляет собой некоторый динамический эффект, возможные 
механизмы возникновения которого на уровне разных деревьев в настоящее время 
дискутируются с биологами. 

В связи с тем, что синхронный кластер охватывает большую группу удаленных 
структурных элементов почки, дополнительная информация об их динамике может 
быть получена путем анализа мгновенных фаз колебаний. Такой анализ позволяет, 
в частности, оценить фазовый сдвиг между колебаниями нефронов, находящихся на 
расстоянии друг от друга. Как отмечалось в работе [21], для соседних нефронов ти- 
пичен режим синфазной синхронизации, который наблюдается более чем в 90% экс- 
периментов. Любопытно отметить, что и для удаленных нефронов, принадлежащих 
единому синхронному кластеру, характерен эффект синфазной (или почти синфаз- 
ной) динамики. В связи с тем, что представление разностей фаз для больших ансам- 
блей автоколебательных систем является не очень наглядным, на рис. 4 продемон- 
стрированы временные зависимости действительных частей вейвлет-коэффициентов 
(это представление аналогично А соз ф в классическом спектральном анализе). В со- 
ответствии с рисунком, режим синфазной синхронизации наблюдается лишь для от- 
дельных участков экспериментальных записей: несинхронная динамика значитель- 
ной части нефронов сменяется синхронной и наоборот, то есть мы наблюдаем эф- 
фект перемежающейся синхронизации. Для наглядности на рисунке приведены за- 
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Рис. 4. Временные зависимости действительных частей вейвлет-коэффициентов, позволяющие проил- 
люстрировать переход от режима несинхронных колебаний значительной части корковых нефронов (а) 
к режиму синфазной синхронизации большой группы структурных элементов почки (б) 

висимости действительных частей вейвлет-коэффициентов только для тех нефронов, 
которые принимают участие в формировании кластера. Соответствующие зависимо- 
сти для нефронов, не демонстрирующих режим синхронных колебаний, не приво- 
дятся, чтобы не загромождать рисунок. 

3. Заключение 

В рамках данной работы проводилось исследование коллективной динамики 
больших групп нефронов, расположенных на поверхности почки. В ходе экспери- 
ментов, проводившихся в Институте Панум (Университет Копенгагена), были за- 
регистрированы сигналы, отражающие скорость кровотока в артериолах нескольких 
десятков корковых нефронов. Далее с помощью метода вейвлет-анализа осуществля- 
лось выделение мгновенных частот ритмов колебаний, обусловленных механизмом 
КГОС. Построение полученных частотных зависимостей на одном графике позволя- 
ет определить участки времени, в течение которых часть нефронов синхронизирует 
свои колебания, и мгновенные частоты соответствующих ритмов подстраиваются. 
В ходе проведенных исследований было установлено, что в формировании класте- 
ров синхронной динамики корковых нефронов участвуют структурные элементы, 
относящиеся к разным нефронным «деревьям», а захват частот и фаз колебаний в 
коллективной динамике нефронов происходит на определенных участках экспери- 
ментальных записей. При этом было обнаружено, что значительные группы нефро- 
нов, расположенные на разных участках поверхности почки, демонстрируют эффект 
синфазной синхронизации колебаний. 

Важно отметить, что сам по себе эффект формирования кластера синхрони- 
зации большой группы корковых нефронов не является случайным наблюдением - 
этот вывод сделан на основе анализа данных экспериментов на пятнадцати крысах, и 
в каждом случае фиксировалось наличие синхронной динамики по крайней мере для 
десяти процентов нефронов, расположенных на поверхности почки. Остается непо- 
нятным, участвуют ли в формировании синхронного кластера глубинные нефроны? 
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Для них характерны более длинные петли обратной связи, что сопровождается более 
длительным движением фильтрата по канальцам и меньшими значениями частоты 
колебаний, связанных с механизмом канальцево-гломерулярной обратной связи. Не 
исключено, что глубинные нефроны демонстрируют режимы синхронной динамики 
с несколько отличающимися периодами. Возможным косвенным подтверждением 
этого предположения служит макроскопическая динамика почки, например, записи 
величины потока крови в почечных артериях. В данных сигналах можно идентифи- 
цировать отдельные ритмические составляющие (по крайней мере при локальном 
спектральном анализе с применением вейвлет-преобразования) [35], но если бы и 
корковые, и глубинные нефроны почки принимали участие в формирование едино- 
го кластера, следовало бы ожидать более четкого проявления характерного пика в 
спектре сигнала, отражающего макродинамику почки. 

Авторы выражают искреннюю признательность N.-13. I ІоІзІеіп-РаіЫои и 
ЭЛ. МагзЬ за предоставленные данные и многочисленные дискуссии. 

Проводимые исследования были поддержаны Министерством образования и 
науки РФ в рамках ФЦП «Научные и научно-педагогические кадры инновационной 
России» на 2009-2013 годы (госконтракт П451). 
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8УМСШЮМ2АТЮ1Ѵ ОГ 08СІІХАТКЖ8 
Ш ТНЕ БУ1ЧАМІС8 ОР Е1Ч8ЕМВЕЕ8 ОГ 8 11 К РА СЕ 1ЧЕРШЮ1Ч8 

О. N. Раѵіоѵа, АЖ. Раѵіоѵ, А.А. Апштоѵ, 

А.І. Магітоѵ, О. V. 8отоѵиеѵа 

Вазей оп Йіе апа1у8І8 о Г ехрегітепіаі йаіа ѵѵс 8Іийу Йіе сойесііѵе йупатісз о Г 
еп8етЫе8 Ггот зеѵегаі 1еп8 перЬгопз Іосаіей оп а кійпеу зшТасе. Изіп§ ѵѵаѵсісі-апаіузіз, 
ІЬе рЬепотепоп о Г 1оскіп§ оГ іп8Іап1апеои8 Ггсуиспсісз апй рЬазез І8 зіийісй іЬа! І8 
саизей Ьу Йіе 1иЬи1о-§1отет1аг ГеейЬаск. И І8 зЬоѵѵп 1Ьа1 зігисіигаі ипіІ8 оГ Йіе кійпеу 
геіаіей Іо йізііпсі перЬгоп 1гее8 рагіісіраіе іп сіизіегз Гогпіаііоп. ТЬе епігаіптепі о Г 
Ггсуиспсісз апй рЬазез оГ 08СІ11айоп8 Гог 1аг§е §гоирз оГ перЬгопз оссигз опіу Гог 8оте 
Гга§тепіз оГ ехрегітепіаі йаіа. И І8 8Іа1ей 1Ьа1 зі§пійсапі §гоирз оГ перЬгопз ріасей іп 
йіГГегеп! агеаз оГкійпеу зигГасе йетопзігаіе ІЬе рЬепотепоп оГіп-рЬазе зупсЬгопігаііоп. 

Кеу\ѵопѣ: ѴсрЬгопз, гепаі Ыоой-йоіѵ аи1оге§и1айоп, ѵѵаѵсісі-апаіузіз, зупсЬгопігаііоп. 


23 






Павлова Ольга Николаевна - родилась в Саратове (1985). Окончила физи- 
ческий факультет Саратовского университета по специальности биохимическая 
физика (2007). В настоящее время работает ассистентом на кафедре радиофи- 
зики и нелинейной динамики СГУ. Область научных интересов: динамика био- 
логических систем. Имеет около 20 работ, опубликованных в отечественной и 
зарубежной печати. 

410012 Саратов, ул. Астраханская, 83 

Саратовский государственный университет им. Н.Г. Чернышевского 
Е-таі1: раѵ1оѵа_о1уа@уайоо.сот 


Анисимов Алексей Александрович - аспирант кафедры радиофизики и нели- 
нейной динамики физического факультета СГУ (специальность «биофизика»). 
Инженер кафедры. Область научных интересов: динамика биологических си- 
стем, анализ временных рядов. 

410012 Саратов, ул. Астраханская, 83 

Саратовский государственный университет им. Н.Г. Чернышевского 
Е-таі1: а1ехеу.а.апІ8Ітоѵ@8іпаі1.сот 


Назимов Алексей Игоревич - студент 5-го курса физического факультета 
СГУ. Область научных интересов: методы анализа структуры сигналов, дина- 
мика биологических систем. Соавтор 3 статей по данной тематике. 

410012 Саратов, ул. Астраханская, 83 

Саратовский государственный университет им. Н.Г. Чернышевского 
Е-таі1: Ка2Ітоѵ191288@уашіех.ги 



Сосновцева Ольга Владимировна родилась в 1966 году. Окончила физи- 
ческий факультет Саратовского университета (1989). Защитила кандидатскую 
диссертацию (1996). Неоднократно работала за рубежом: в Германии, Дании и 
Южной Корее. В настоящее время работает в Датском техническом универси- 
тете (г. Люнгбю). Научные интересы - явления самоорганизации в ансамблях 
нелинейных хаотических и стохастических систем; особенности синхрониза- 
ции хаотической динамики и стохастических процессов, динамика биологиче- 
ских систем. Имеет около 70 публикаций в отечественной и зарубежной печати. 

Оептагк 1017 СорепЬа§еп, №>гге§а1е 10, Ро8іЪох 2177 
Рагшт ІпзІіШіе, Ііпіѵегзііу оГ СорепЬа§еп 
Е-таі1: 0І§а@8ітс1.ки.с1к 


24 




